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ABSTRACT 

We find the orbit of the Neptune-sized exoplanet HAT-P-1 lb to be highly inclined relative to the equatorial 
plane of its host star. This conclusion is based on spectroscopic observations of two transits, which allowed 
the Rossiter-McLaughlin effect to be detected with an amplitude of 1.5 m s" 1 . The sky -projected obliquity 
is 103!^ de grees. This is the smallest exoplanet for which spin-orbit alignment has been measured. The 
result favors a migration scenario involving few-body interactions followed by tidal dissipation. This finding 
also conforms with the pattern that the systems with the weakest tidal interactions have the widest spread in 
obliquities. We predict that the high obliquity of HAT-P-11 will be manifest in transit light curves from the 
Kepler spacecraft: starspot-crossing anomalies will recur at most once per stellar rotation period, rather than 
once per orbital period as they would for a well-aligned system. 

Subject headings: planetary systems — planets and satellites: formation — planet-star interactions — stars: 
rotation 



1, INTRODUCTION 

The origin of close-in planets is the longest-standing prob- 
lem in exoplanetary science (Mayor & Queloz 1995). Re- 
cently, the orbits of some close-in planets were found to be 
highly inclined relative to the equatorial planes of their host 
stars (see, e.g., Hebrard et al. 2008, Narita et al. 2009, Winn et 
al. 2009, Triaud et al. 2010). This evidence supports theories 
for close-in planets in which their orbits shrink due to gravi- 
tational perturbations from other bodies followed by tidal dis- 
sipation (Matsumura et al. 2010). The evidence disfavors the 
other leading theory, in which the orbits shrink due to grad- 
ual interactions with the protoplanetary gas disk, unless the 
disks were somehow misaligned with their host stars (Bate et 
al. 2010, Lai et al. 2010). 

To this point, spin-orbit alignment has been measured only 
for "hot Jupiters," with masses ranging from 0.4-20 Mj up . We 
would like to extend these studies to smaller planets, in order 
to see whether they migrate in a similar way as larger plan- 
ets, and to understand which factors are associated with or- 
bital misalignment. It has been claimed, for example, that 
tilted orbits are more prevalent for massive planets (Johnson 
et al. 2009, Hebrard et al. 2010), or for stars with thinner con- 
vection zones (Winn et al. 2010, Schlaufman 2010). 

Here we present a spin-orbit study of HAT-P-1 lb, a "hot 
Neptune" of mass 0.08 Mj up and radius 0.42 Rj up on an ec- 
centric, 4.9-day orbit around a K4 dwarf (Bakos et al. 2010; 
B10 hereafter). We observed the Rossiter-McLaughlin (RM) 
effect (§ 2), and modeled it (§ 3), finding the orbit and stellar 
spin to be misaligned (§ 4). 
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2. OBSERVATIONS AND DATA REDUCTION 

We obtained 132 new spectra of HAT-P-11 with the High 
Resolution Spectrograph (HIRES; Vogt et al. 1994) on the 
Keck 1 10m telescope. Most of the new spectra were gathered 
on nights when transits were predicted. On 2009 Aug 2/3 
we gathered 7 spectra during a transit, although fog pre- 
vented us from observing before or after the transit. On 
2010 May 26/27 we obtained 32 spectra starting at around 
first contact and extending for a few hours beyond the tran- 
sit. On 2010 Aug 22/23 we obtained 70 spectra spanning the 
entire transit and a few hours beforehand and afterward. The 
remaining 23 spectra were obtained sporadically throughout 
the 2009-2010 observing season. 

We used the instrument settings and observing procedures 
that are standard for the California Planet Search (Howard et 
al. 2009). In particular, we used an iodine gas absorption cell 
to track the instrumental response and wavelength scale. The 
radial velocity (RV) of each spectrum was measured with re- 
spect to an iodine-free template spectrum, using a descendant 
of the algorithm of Butler et al. (2006). Measurement errors 
were estimated from the scatter among the fits to individual 
spectral segments spanning a few Angstroms. Table 1 gives 
all the Keck/HIRES RVs, including re-reductions of the 50 
spectra presented by B10. 



3. analysis 

Merging all the RVs into a single analysis requires some 
care because the host star is chromospherically active. B10 
found a photometric signal with period 29.2 days and ampli- 
tude 3 mmag, which they attributed to starspots being carried 
around by stellar rotation. One would expect a corresponding 
RV signal at the same period and its harmonics, with an am- 
plitude of order 0.3% of the projected stellar rotation speed 
(vsin/i), or approximately 5 ms" 1 . Indeed, B10 found evi- 
dence for "stellar jitter" of amplitude 5 ms -1 , supporting this 
interpretation. 

We investigated this issue by fitting the out-of-transit RVs 
with a model consisting of a single Keplerian orbit plus a con- 



2 



Winn et al. 2010 



stant acceleration^ and seeking evidence for time-correlated 
residuals. As seen in Figure 1, the residuals are strongly cor- 
related on timescales shorter than 5-10 days, as expected. On 
longer timescales there are no obvious correlations. 
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Fig. 1. — Correlations of the RV residuals. Top. — Products of pairs of 
residuals, as a function of the time elapsed between the observations. Sig- 
nificant positive correlations are seen for At ^ 5 days. Bottom. — Same, but 
restricting the analysis to only one data point per observing run, giving a 
minimum time separation of 20 days. No significant correlations are seen. 

While it would be possible to model the RV covariances, we 
chose the simpler approach of selecting a subset of RVs that 
are effectively independent. Specifically, we chose a single 
spectrum from each observing run, resulting in a sample of 
21 out-of-transit RVs spaced apart by a minimum of 20 days. 

We turn now to the transit nights. During a single night, 
one would expect the rotationally modulated RV signal to act 
as a nearly constant offset. Therefore, each transit night was 
assigned a free parameter that shifts all the RVs by a common 
amount; only the intranight variations were deemed signifi- 
cant. With this approach the data from 2009 Aug 2/3 were 
rendered useless, because no data were gathered outside of 
the transit. We omitted those data from further consideration. 

Our model for the RVs was the sum of the Keplerian orbital 
motion, a constant acceleration, the RM effect, and the offsets 
described above. The fitting statistic was 
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where the first term is the usual sum of squared residuals, and 
the other terms represent a priori constraints on parameters 
that were determined more precisely from the larger body of 



data analyzed by B 10. In this expression Pjays is the orbital pe- 
riod in days, BJD C is a particular time of inferior conjunction; 
Tdays is the time between first and fourth contact; Td ays is the 
time between first and second contact; R p and are the radii 
of the planet and star; and vsim* is the star's sky -projected 
rotation speed. 

Each of the 21 orbital RVs was assigned an error bar ct, 
equal to the quadrature sum of the measurement error and a 
"jitter" of 5.5 m s , the value giving \ 2 = Ndof when the or- 
bital RVs were fitted alone. For the transit-night RVs, the jitter 
was fixed by the requirement \ 2 = AW when fitting the data 
from that night along with the orbital RVs. The results were 
1.8 and 1.5 m s" 1 for 2010 May 26/27 and 2010 Aug 22/23, 
respectively. The relative smallness of these values corrobo- 
rates our assumption that the activity-induced RV variations 
occur mainly on longer timescales. A similar contrast be- 
tween intranight and internight jitter was observed previously 
for HD 189733, another active K star (Winn et al. 2006). 

We modeled the RM effect with the technique described by 
Winn et al. (2005), finding in this case that a sufficiently ac- 
curate description for the RV shift is the product of the loss of 
light and the RV of the portion of the stellar photosphere be- 
neath the planet. We neglected differential rotation, and took 
the stellar limb-darkening law to be linear with a coefficient 
of 0.79 (Claret 2004). 

Parameter optimization and error estimation were achieved 
with a Markov Chain Monte Carlo algorithm, using Gibbs 
sampling and Metropolis-Hastings stepping. Table 2 gives 
the results for each parameter, based on the 15.85%, 50%, and 
84. 15% confidence levels of the marginalized a posteriori dis- 
tributions. Figure 2 shows the RV data: the left panel shows 
the orbital RVs; and the right panel shows the transit-night 
RVs after subtracting the calculated variation due to orbital 
motion, thereby isolating the "anomalous RV" due to the RM 
effect. Figure 3 shows the a posteriori distributions for the 
key parameters vsin/* and A (the sky position angle from the 
stellar north pole to the orbital north pole). 

The most important result is A = 103^q degrees, indicating 
a major misalignment between the stellar rotation axis and the 
orbit normal. Qualitatively this follows from the observation 
that the RM effect was observed to be a blueshift throughout 
the transit, as opposed to the "red-then-blue" pattern of a well- 
aligned system. 

4. DISCUSSION 

Because the signal has an amplitude of only 1.5 m s" 1 , 
smaller than any other RM signal yet reported, it is impor- 
tant to test the robustness of the results. First, we tried an- 
alyzing each transit individually, rather than combining the 
data from both transits. As shown in Figure 3, the results 
are in good agreement^ Second, we repeated the analysis 
without the a priori constraint on vsin/ + , the other parame- 
ter of greatest relevance to the RM effect. The results were 
vsin/ + = 1.13^07Q km s" 1 and A = 100^| 8 degrees. Third, we 
checked for any correlations between the RM signal and the 
strength of Ca II H and K emission, or the shape parameters 
of the instrumental line spread function. Significant corre- 
lations would have raised suspicion of systematic errors, but 
none were found. 



7 The evidence for a constant acceleration (linear velocity trend), and the 
implied existence of another orbiting body besides HAT-P-llb, were estab- 
lished by B10. 



8 In fact there are three mutually reinforcing datasets: Hirano et al. have 
submitted a paper reporting A = 103+^ degrees, based on independent obser- 
vations of HAT-P-llb (T. Hirano and N. Narita, private communication). 
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FIG. 2. — Radial velocities and the RM effect. Left. — Spectroscopic orbit of HAT-P-1 1, based on the subset of 21 RVs analyzed here. Right. — Anomalous 
RV of HAT-P-1 1 spanning two transits (top and bottom), and the time-binned combination (middle; binned x7 with a maximum bin size of 0.5 hr). The orbital 
contribution to the RV has been subtracted. The solid line shows the best-fitting model of the RM effect. The dashed curve shows the best-fitting "well-aligned" 
model (A = 0, vsini* = 1.3 km s -1 ), which is ruled out with A\ 2 = 52.4. 



In addition, if the 29.2 day periodicity detected by B10 is 
indeed the rotation period, then a powerful test for spin-orbit 
misalignment is available. If the system were well-aligned we 
would have A = and z* w 90°. This would imply 



271-/?* 



vsim* = v = 



= 1.30 km s" 



(2) 



where we have used = 0.752 Rq and P rot = 29.2 days (B10). 
When refitting the data with these constraints, the minimum 
X 2 rises from 1 1 1.9 to 164.3 (Ax 2 = 52.4), with 1 14 degrees 
of freedom. Thus the well-aligned model is ruled out with 
7.2<7 confidence: either A is large, or else z* must be far from 
90° to be compatible with the low amplitude of the RM signal. 
The best-fitting well-aligned model is illustrated with a gray 
dashed curve in Figure 3. 

HAT-P-1 lb is the first "hot Neptune" for which the RM ef- 
fect has been measured. Our results suggest that tilted orbits 
are common for hot Neptunes, just as has been found for hot 
Jupiters. The same migration mechanisms that are invoked 
to explain the larger planets with tilted orbits — gravitational 
scattering by planets, or the three-body Kozai effect — may 
also have operated in this case. It should be noted that the 
spin-orbit results are not the only evidence for a perturba- 
tive origin for many close-in planets. Further evidence comes 
from their occasionally high orbital eccentricities, the clus- 
tering of their orbital distances near the value expected from 
tidal circularization, and their tendency to lack companions 
with periods between 10-100 days (Matsumura et al. 2010). 

Since HAT-P-1 lb is the lowest-mass planet yet probed by 
RM measurements, and it is misaligned, our findings are at 
odds with the hypothesis that misalignments occur mainly for 



the most massive planets (Johnson et al. 2009). They do, how- 
ever, support the correlation between large obliquity and or- 
bital eccentricity (Johnson et al. 2009, Hebrard et al. 2010), 
as the orbit of HAT-P-1 lb has a significant eccentricity. 

Another emerging trend is that misalignments occur mainly 
for stars with high effective temperatures or large masses 
(7 eff > 6250 K or M* > 1.2 M Q ). Winn et al. (2010) spec- 
ulated that this is due to tidal interactions: cool stars realign 
with the orbits, but hot stars cannot realign because tidal dissi- 
pation is weaker in their thinner outer convection zones. The 
HAT-P-11 system is an important test case because the star 
is cool and low-mass, and yet tidal interactions are weak due 
to the planet's relatively small size and long period. If stellar 
temperature or mass are the determinants then one would ex- 
pect HAT-P-1 1 to be well-aligned like other cool stars. But if 
tides are the underlying factor, then HAT-P-1 1 would be mis- 
aligned, as we have observed. Specifically, with reference to 
Eqn. (2) of Winn et al. (2010), HAT-P-11 experiences even 
weaker tides than WASP-8, a cool star already known to have 
a high obliquity0 Thus, HAT-P-11 is a telling exception to 
the rule that hot stars have high obliquities: it implicates tidal 
evolution as the reason for low obliquities among cool stars 
with more massive planets in tighter orbits. 

By good fortune, HAT-P-11 is in the field of view of the 
Kepler spacecraft (Borucki et al. 2010). The precise photo- 
metric time series will allow the candidate 29.2-day rotation 
period to be checked. Asteroseismological studies may reveal 
the stellar mean density, age, inclination, and other parame- 

9 Matsumura et al. (2010) have also argued that tidal dissipation in the 
HAT-P- 1 1 system would be too weak to realign the star with the orbit, despite 
the star's thick outer convective zone. 
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FIG. 3. — Results for the RM parameters. Left. — In this polar plot, the angular coordinate is A and the radial coordinate is p(X), the marginalized posterior 
probability distribution. Results are shown for analyses based only on individual transits, as well as for the entire dataset. Right. — Joint constraints on A and 
vsirn*. The contours represent 68%, 95%, and 99.73% confidence limits. The marginalized posterior probability distributions are shown on the sides of the 
contour plot. The dashed histogram shows the a priori constraint on v sin U . 



ters (Christensen-Dalsgaard et al. 2010). Furthermore we pre- 
dict that Kepler will see a pattern of anomalies in the transit 
light curves that will betray the system's spin-orbit misalign- 
ment. As usual for a spotted star, there will be a "bump" or 
"rebrightening" in the transit light curve whenever the planet 
occults a starspot (see, e.g., Rabus et al. 2009). For a well- 
aligned star, the bumps recur in successive transits for as long 
as the spot is on the visible hemisphere. There is a steady ad- 
vance in phase of the bumps due to the star's rotation between 
transits. However, for a star like HAT-P-1 1 with A ss 90°, the 
events will not recur in this manner, because the star's rotation 
moves the spot away from the transit chord. A spot must com- 
plete a full rotation before returning to the transit chord, and 
even then, the planet will miss it unless it has also completed 
an integral number of orbits. 

For HAT-P-1 1, it happens that P ro t/P ib ~ 6. If the star were 
well-aligned and had one spot initially on the transit chord, 
then we would typically see an alternation between 2-3 light 
curves with bumps, and 2-3 without bumps (when the spot 
is on the far side). But because of the misalignment, spot 



anomalies will only recur every 29.2 days, after the star has 
rotated once and the planet has completed 6 orbits. Complica- 
tions may arise due to differential rotation, as well as the mul- 
tiplicity and evolution of spot patterns. Nevertheless, this phe- 
nomenon should allow for an independent test of spin-orbit 
misalignment for HAT-P-1 1 as well as other spotted stars. 
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TABLE 1 

Relative Radial Velocity Measurements of HAT-P- 1 1 



BJD UTC RV[ms _1 ] Error [ms _1 ] 

2454335.89332 -1.33 1.01 

2454335.89998 -2.66 1.03 

2454336.74876 -1.30 0.90 

2454336.86163 -4.74 1.02 

2454336.94961 -7.56 0.94 

NOTE. — The RV was measured relative to an arbitrary template spectrum; only the differences are significant. The uncertainty given in Column 3 is the internal error only and does 
not account for any possible "stellar jitter." (We intend for this table to be available in its entirety in a machine-readable form in the online journal. A portion is shown here for guidance 
regarding its form and content.) 



TABLE 2 
Model Parameters 

Parameter Value 
Parameters controlled mainly by priors 

Orbital period, P [days] 4.88781501 ± 0.0000068 

Midtransit time [BJD UTC ] 2,454,605.89130±0.00032 

Planet-to-star radius ratio, R,,/R* 0.0576 ± 0.00090 

Orbital inclination, i [deg] 89.17^ 

Fractional stellar radius, R+/a 0.0673 ± 0.0018 

Parameters controlled mainly by RV data 

Velocity semiamplitude, K [m s _1 ] 12.9 ±1.4 

ecosw 0.261 ±0.082 

esinw 0.085 ±0.043 

RV offset, 2010 May 26/27 [m s _1 ] -19.8 ± 3.9 

RV offset, 2010 Aug 22/23 [m s" 1 ] -17.5 ±4.2 

RV offset, all other data [m s -1 ] -13.0 ± 2.2 

RV trend, 7 [m s~' day* 1 ] 0.0185 ± 0.0036 



Projected stellar rotation rate, vsinf* [km s '] 1.00; 
Projected spin-orbit angle, A [deg] 103 



+0.95 
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